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Abstract 



^ r~| A superposition rule for two solutions of a Milne-Pinney equation 



is derived. 



1 Introduction 

r- -n ■ 
>, 

The search for explicit solutions of systems of first-order differential equa- 
tions is not an easy task and there exist only few cases in which we can 
express in a compact way the form of the general solution. The determina- 
te ■ tion of particular solutions can be very useful in the reduction of the original 
problem to a simpler one. In the simplest cases of linear systems the linear 
Q ! superposition principle plays a relevant role and we are able to express the 
general solution as an arbitrary linear combination with real coefficients of a 
fundamental set of solutions. There are also cases of nonlinear (systems of) 
differential equations for which a nonlinear superposition rule does exist, the 
best known one being the Riccati equation 

dx 

— = b (t)+b 1 (t)x + b 2 (t)x\ (1) 

which is very often used in many fields of mathematics, control theory and 
theoretical physics (see for instance [H [2] and references therein) and its 
importance in this field has been increasing since Witten's introduction of 
super symmetric Quantum Mechanics [3]. 



In fact, if three particular solutions of ([T]) are known, x±,x 2 and x 3 , the 
general solution can be found from the cross ratio relation 

x - x 2 x 3 - x 2 
which provides us a nonlinear superposition rule: 

kx±(x 3 - x 2 ) + x 2 {xi - x 3 ) 



X = (/)(x 1 ,X2,X 3 ,k) 



k (x 3 - x 2 ) + (xi - x 3 ) 



The general theory for nonlinear superposition rules was developed by 
Lie [I] and it has been more recently revisited from a geometric perspective 

W 

Another nonlinear equation with relevance in physics is the today called 
Milne-Pinney equation: 

k 
x = -uj 2 {t)x + — , A;>0. (2) 



.r 



This equation is defined on R — {0} and it is invariant under parity, i.e. if 
x(t) is a solution then x(t) = —x(t) is a solution too. That means that it is 
enough to restrict ourselves to the x > case. 

Equation ([2]) models many physical problems such as propagation of laser 
beams in nonlinear media, plasma dynamics, or the mean field dynamics for 
Bose-Einstein condensates through the so-called Gross-Pitaevskii equation. 
It was introduced by Ermakov as a way for finding a first-integral of the 
corresponding time-dependent harmonic oscillator described by 

y + u 2 (t)y = (3) 

by means of the Ermakov system [7] made up by the two equations ([2]) and 
((3]). This Ermakov system admits the well-known Ermakov invariant 

(7/ \ 2 
-J + (xv y - yv x ) 2 , (4) 

from where one can see that the general solution x(t) of Milne-Pinney equa- 
tion can be expressed in terms of two solutions y\ and y 2 of the associated 
harmonic oscillator as follows [8]: 

x = ^-Jcwf + C 2 y\ ± fe - kW* yi y 2 , (5) 



where W is the Wronskian of the two solutions of fl3]), W = 2/1^2 — 2/22/1, which 
is constant, and C\ and C 2 are non-negative constants such that \C\C 2 > 
kW 2 . In this way, as 

(Vc~m - Vc~2V2) 2 > 



then Ciyf + C 2 y\ > V^C^C^y 1 y 2 > \I^C\C 2 - kW 2 y 1 y 2 and thus the solu- 
tions given by flS} are real. 

In a similar way, it has recently been shown in [U] that there exists also a 
superposition rule involving three solutions of a related Riccati equation for 
giving the general solution of (J2J). 

A generalisation of this Ermakov system was proposed in PJJ [TTJ [121 EE] : 



x = -^f{y/x)-uj 2 {t)x 



x c 



(6) 



y = -^g{y/x) - u 2 {t)y 

and a new invariant was obtained for this generalised Ermakov system: 



I(x, y, v x , v y ) = (xv y - yv x ) 2 + 2 
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KuJ_ 



du . (7) 



This generalised Ermakov system reduces to the standard one for f(u) = k 
and g(u) = 0. In this case ((7|) becomes (j4j). 



2 A new superposition rule for the Milne— 
Pinney equation. 

Our aim in this section is to show that there exists a superposition rule for 
the Milne-Pinney equation (J2J) [HI [HI [IS] in terms of a pair of its particular 
solutions. 

In fact, one see from ((Tj) that in the particular case of / = g = k, if a 
particular solution xi is known, there is a ^-dependent constant of motion for 
the Milne-Pinney equation given by (see e.g. [To]): 



I\ = (x\x — x\x) 2 + k 



-) + 

X\ I \ X / 



EiV 



(8) 



If another particular solution x 2 of © is given, then we have another 
^-dependent constant of motion 



(x 2 x — x 2 x) 2 + k 



x 

X2 



X2\' 
X J 



(9) 



Moreover, the two solutions x\ and x 2 provide us a function of t which is 
constant and generalises the Wronskian W of two solutions of ([2]): 



(#1X2 - x 2 xi) 2 + k 



X2 
Xi 



■1-2 



(10) 



Remark that for any real number a the inequality (ft — 1/ft) 2 > implies 

o 1 



a 



> 2. 



a" 



and the equality sign is valid if and only if \a\ 



ft 2 4 



a 



1 



a' 



Therefore, as we have considered k > 0, we see that Jj > 2 k, for i = 1, 2, 3. 
Moreover, as the solutions x\(t) and x 2 (t) are different solutions of the Milne- 
Pinney equation, it turns out that ^3 > 2k. 

The knowledge of the two constants of motion Ii and I 2 , together with 
the constant value of ^3 for a pair of solutions of ([2]), can be used to obtain 
the superposition rule for such a differential equation. In fact, given two 
solutions x\ and x 2 of (T5]), first the constant of motion Qj allows us to write 
an explicit expression for x in terms of x, x 2 and I 2 . 



x 




k— + / 2 -2 - k^ 



,t n 



,! .■ ;- 



X 



2 • 



and using such expression in ([8]) we see after a careful computation that x 2 
satisfies the following fourth degree equation 



(If - 4fc 2 



).?'-, 



2(/ 1 / 2 - 2I 3 k)xix 2 2 + (It 



Ak 2 )xi- 



- 2{{hh - 2hk)x\ + (J1J3 - 2I 2 k)x\)x 2 + (/ 3 2 - 4fc 2 )x 4 = , (11 



where use has been made of the constant value of / 3 for two solutions of the 
Milne-Pinney equation. 

Therefore, we can obtain from (fTTj) the expression for the square of the 
solutions of ([2]) in terms of any pair of its particular positive solutions by 
means of a superposition rule 



x 2 = \ lX \ + \ 2 x\ ± 2^\ 12 [-k(xf + xj) + hx\x 2 2 l (12) 

where the constants Ai and A2 are given by 

Al " I 2 -4k 2 ' X2 ~ I 2 -4k 2 ' (13) 

and A12 is a constant which depends on Ai, A2, k and the constant value of 
7 3 as follows: 

A 12 = A 12 (A l5 A 2 ; I 3 , k) = Ws + M-l^+Aj + Al) = ^ /2; ^ k) (14) 
where the function (p is given by 

„ {I J-J / ,_ /l/2/3-(/ 1 2 + / 2 2 + / 3 2 )fc + 4fc3 

y^i, V 2 , V 3 , fcj - 2 - 2 . (15) 

It is important to remark that if Ai < then A 2 > and if A 2 < then 
Ai > 0, i.e. if Ai < then J 2 / 3 < 2Iik, and thus J 2 < 2kl\j I3. Therefore, 
A 2 (J 3 2 - 4k 2 ) = hh ~ 2kl 2 > hh - 4k 2 h/I 3 = h(I 2 - 4k 2 ) > 0, and thus, 
as I 3 > 2k, A 2 > 0. Similarly we obtain that A 2 < implies Ai > 0. 

The parity invariance of ([2]) is displayed by (fl2"l) which gives us the solu- 
tions 

x 2 = Aix 2 + X 2 x 2 2 ± 2^/\ l2 [-k(xj + xj) + I 3 x\x\ ] . (16) 

In order for x 2 to be a positive real number giving rise to a real solution 
of the Milne-Pinney equation the constants Ai and A 2 in the preceding ex- 
pression should satisfy some additional restrictions. In particular, they must 
be such that 

Ai2[-fc(zi(0) + 4(0)) + / 3 * 2 (0)x 2 (0) ] > 
and 



Aix 2 (0) + A 2 x 2 (0) ± 2^/ \ 12 [-k(xf(0) + xl(0)) + l 3 x 2 (0)x 2 (0)} > 0. 



If these conditions are satisfied, then differentiating (1161) in t = for x\ = 
Xi(t) and x 2 = x 2 (t) solutions of fl2]), it can be checked out that x(0) is also 
a real constant. As x(t) is a solution with real initial conditions then x(t) 
given by ([TBI) is real in an interval of t and thus all the obtained conditions 
will be valid in an interval of t. 

If we take into account that we have consider x 2 > 0, we can simplify the 
study of such restrictions by writing (TIT)]) in terms of the variables x 2 and 
z = (xi/x 2 ) 2 as 



.2 



Ju Juty 



[\ x z + A 2 ± 2v/Ai 2 [-fc(z 2 + l) + / 3 z]) , (17) 

and the preceding conditions turn out to be Ai 2 [— k(z 2 + 1) + I 3 z] > and 
X lZ + A 2 ± 2 v /\ 12 [-k(z 2 + l) + I 3 z] > 0. 

Next, in order to get X 12 [—k(z 2 + 1) + I 3 z] > 0, we first note that this 
expression is not definite because its discriminant is A^ 2 (/f — 4/c 2 ) > 0, and 
this restricts the possible values of Ai and A 2 for a given z. With this aim 
we define the polynomial P(z) given by 

P(z) = -k(z 2 + l) + I 3Zj (18) 

with roots 

, = ^ = /,*^Fg (19) 

which can be written in terms of the variable a 3 = I 3 /2k as 



z± = a 3 ±Ja 2 3 -l. (20) 



As a 3 > 1, then a 3 > \Ja 3 — 1 > and thus z± > 0. The sign of the 
polynomial P(z) is displayed in Fig. 1. 

The region M.+ x M + splits into three regions, 

A = {(x 1 ,x 2 ) E l + xl + | x\ > t/z^x 2 } [J{(xi,x 2 ) G M + xlR + | x x < y/zZx 2 } 



B = {(xi,x 2 ) G M + x M + | ^pz~Zx 2 < x\ < v /ZjTx 2 } 
separated by the region 

C = {(xi,x 2 ) G M + xlR + | x\ = y/z^x 2 ] |){(xi,x 2 ) G 1R + xM + | x\ = v / 5Tx 2 } 



Xl 




«3^2 



r = Xi 



Z-X 2 



■1*2 



Figure 1: Sign of the polynomial P(xi,X2). 



of the straight lines X\ = JY^Xi and x\ = y/zlx2- The condition to make 
Ai2-P(z) non-negative in region A where P takes negative values is to choose 
Ai and A2 such that Ai2(Ai, A2, ^3, k) < 0. Similarly, as P is positive in region 
B we have to choose Ai and A2 such that Ai2(Ai, A2, ^3, k) > 0. Finally as P 
vanishes in region C there is no restriction on the coefficients Ai and A2. 

Once we have stated the conditions for Ai2-P(z) to be non-negative we 
still have to impose the condition 



X x z + A 2 ± 2^X 12 [-k(z 2 + l) + hz] > 0. 



(21) 



In order to study these conditions we study the sign of the polynomial 
P hyk {z, Ai, A 2 ) = (Aiz + A 2 ) 2 - AX 12 [-k(z 2 + 1) + hz\ 



where 



+ (aXi + b\ 2 ) 



/f - Ak 2 



AP(z)k , 



I 2 - Ak 2 



+ z z 



AP(z)k 
I 2 - Ak 2 ' 



(22) 



(23) 



As we remarked before, Ai,A2 cannot be both negative. Let K denote 
the set K = R 2 — {(Ai, A 2 ) G R 2 | Ai < 0, A 2 < 0} and consider three cases: 

1. If (xi, x 2 ) G A, then as P(z) < 0, it must be Ai 2 < in order to satisfy 
Ai 2 P(z) > 0. In this case, if K\ and K 2 are the sets 



K 1 = l{\ 1 ,\ 2 )eK; 



Ko 



(Ai,A 2 )6A r ; 



4P(^)/ 3 
7|-4A; 2 



> |aAi + 6A 2 | 



7-2 
J 3 



4A; 2 



< laAi 



6A 5 



(24) 



We find the following particular subcases 

(a) If (Ax, A 2 ) G K u then P h , k {z, X u A 2 ) > 0. 

(b) If (Ai, A 2 ) G K 2 then P Jaifc («, Ai, A 2 ) < 0, 

that can be summarised by means of Figure 2. 

A 2 




Figure 2: Sign of the polynomial Pi 3i h(z, Ai, A 2 ) in K. 



2. If (xi,x 2 ) G B, as P(z) is positive, then A12 must also be positive, 
Ai 2 > 0. Thus for (Ai, A 2 ) e #1 U if 2 , P/ Slfc («, Ai, A 2 ) > 0. 

3. If (xi, x 2 ) E C, then for (A l5 A 2 ) e K x U # 2 , P/ 3 , fe (z, Ai, A 2 ) > 0. 

In those cases in which Pj 3t j,(z, X\, A 2 ) > we can assert that \X\Z + A 2 | > 
2\/Ai 2 [— k(z 2 + 1) + I 3 ~z], but we still have to impose that X\Z + A 2 > 
for (12~T1) to be positive. Nevertheless, this is very simple, because if such a 
pair (Ai,A 2 ) does not satisfy X±z + A 2 > 0, the pair of opposite elements 
(— Ai, — A 2 ) actually does it, the other conditions being invariant under the 
change Ai — > —Xi with i = 1,2. 

In those cases in which Pj 3! k(z, Ai, A 2 ) < we can assert that \XiZ + A 2 | < 
2\/Ai 2 [— k(x\ + x 2 ) + I3X1X2 ], and in this case the unique valid superposi- 
tion rule is 

x = \x 2 \ (\iz + A 2 + 2y/X 12 [-k(z 2 + 1) + h A) , (25) 

which is equivalent to 

1/2 



x= I X 1 xi + X 2 x 2 2 + 2^X 12 [-k(xf + xi) + I 3 x 2 2 xj}\ . (26) 

It is important to notice that the given conditions for Ai, A 2 are used to 
obtain real solutions of the Milne-Pinney equation and if no restriction were 
considered we would had obtained real and imaginary solutions for the real 
Milne-Pinney equation. 

Expression f|T6|) provides us a (as far as we know, new) superposition 
rule for the positive solutions of the Pinney equation ([2]) in terms of two of 
its independent particular positive solutions. Therefore, once two particular 
solutions of the equation (J2J) are known, we can write its general solution. 
Note also that, because of the parity symmetry of (J2J), the superposition (TTB1) 
can be used with positive and negative solutions. In all these ways we obtain 
if k > non- vanishing solutions of ([2]). It is also to be remarked that the 
same procedure as before can be repeated in the case k < to obtain similar 
results. 

A similar superposition rule works for negative solutions of Milne-Pinney 
equation (J2J): 

x = - \X lX \ + A 2 x^ ± 2^X l2 {-k{x\ + xj) + J 3 x\x\) J . (27) 
where once again x\ and x 2 are arbitrary solutions. 



3 Relations between the new and the known 
superposition rule. 

We can now compare both superpositions rules (jSJ) and (fT6l) and check that 
actually the latter reduces to the former one when x± and x 2 are obtained 
from solutions y\ and y 2 of the associated harmonic oscillator equation. 

Let 7/1 and y 2 be two solutions of (JSJ) and W its Wronskian. Consider 
the two particular positive solutions of the Milne-Pinney-equation x\ (t) and 
x 2 (t) given by 

Xi 

(28) 
x 2 



V2 

\w\ 
\w\ 


-^C iy 2 + C 2 y 2 , 


-y/c 2 yf + C iy l 


'nally 


impose 


4CiC 2 = kW 2 . 



(29) 

-Z3 given by (ITU]) for the two particular solutions of the Milne-Pinney 
equation can then be expressed as a function of the solutions y\ and y 2 of 
the time-dependent harmonic oscillator and its Wronskian W. After a long 
computation ^3 turns out to be 

_ A(C 2 + op 

J 3 - W2 > ( 30 ) 



and then using the explicit form (128]) of the particular solutions and taking 
into account (130]) in f|T6|) we obtain that 



Aix? + \ 2 x\ ± 2^/\ 12 (-k(xf + ar|) + hx\x 2 2 ) = — (C1A1 + C2A2M 

+ (dA 2 + C 2 Ai)2/ 2 2 ) ± ^^(CiAx + C 2 A 2 )(CiA 2 + C 2 Ai) - fc^ 2 y iy2 . 

(31) 

Consequently, from the superposition rule (TT6|) we recover expression (JSJ): 

/2 I ^ 

x = ——\Jnxyl + /i 2 y| ± y/Afufiz - kW 2 y 1 y 2 , (32) 

10 



where now 

'/i 1 = (CiAi + C 2 A 2 ), 

/i 2 = (CiAa + CsAi). 

Once we have stated the superposition rule we still have to analyse the possi- 
ble values of Ai and A2 that we can use in this case. If we use the expression 
(!30|) we obtain after a short calculation the following values z± 

4C| AC 2 

z+ ~kW^ z ~~mn- (33) 

Now if we write y\ and y\ in terms of x\, x\ and W from the system ( |28l) we 
obtain 

Cl^Cl\-C 2 C 1 2 ){ X xl) = {y\j- (34) 

Therefore, as C 2 > C\ the condition of being y\ and , y 2 being positive is 



\ C 2 x\ > C x x\ 



C x x\ < C 2 x\ 
2 



(35) 



and it is verified if x\jx\ < C 2 /d = ACl/kW 2 = z+ and x\lx\ > d/C 2 = 
4Cf/kW 2 = z_, because of (|29l) . Thus, (xi,x 2 ) E B and therefore the 
only restrictions for Ai, A 2 are A i2 > and \\x\ + A 2 x 2 > 0. Obviously, by 
means of the change of variables (|28j) this last expression is equivalent to 
fiiyf + \i 2 y\ > and thus /ii and /i 2 can not be simultaneously negative. 
Furthermore, Ai 2 (lf — 4k 2 ) = ^ii\\i 2 — kW 2 . As we have said that Ai 2 > 
then 4/iiyU 2 > kW 2 , i.e. [i\\x 2 is positive and thus, /ii and /i 2 are positive. In 
this way we recover the usual constants of the known superposition rule of 
the Milne-Pinney equation in terms of solutions of an harmonic oscillator. 

4 A simple example of the new superposition 
rule. 

As a particularly simple but interesting application, we can use the superpo- 
sition rule ffl6l) for the Milne-Pinney equation in order to obtain the general 
solution of the constant unity frequency case: 

k , ^ 

x = — — x . (36) 

x 6 

11 



As indicated by Pinney in [SJ, the superposition rule (TTBT) gives us the 
general solution in terms of solutions of the corresponding harmonic oscillator 

y = -y- (37) 

Choosing two positive constants C\ and C 2 such that 4CiC 2 = kW 2 and 
C\ < C 2 , we obtain the following two independent particular solutions of the 
Pinney-equation: 

y/2 



J (38) 



If we take, for instance, y\ = cost and y 2 = sint, for which W = 1, the 
following solutions for the Milne-Pinney equation are obtained: 



Xx = v^yCi cos 2 (t) + C 2 sin 2 (t), 
x 2 = V2^C2Cos 2 (t)+dsm 2 (t). 

We can use these two expressions in ( fl6l) for obtaining the general solution 
of the Milne-Pinney equation. If we choose C\ = yk/A and C2 = vk then 
we obtain the general solution for the Milne-Pinney equation with constant 
frequency: 

\fk 
x(t) = ^— (5(Ai + A 2 ) + 3(Ai - A 2 ) cos(2t) 



+ 2v/(4Ai + A 2 )(Ai + 4A 2 ) - 4sin(2t)) 1 / 2 . (39) 

Obviously, as we said before, Ai and A 2 are such that Ai 2 > and \\x\ + 
X 2 x 2 > 0. 



5 Conclusions and Outlook. 

In summary, we have obtained a new superposition rule for the Milne-Pinney 
equation in terms of a set of two of its particular solutions. Also, we have 
related the new superposition rule with the one proposed in [8J for solutions 
of the related equation ([3]). Finally, a particular example in which we use this 

12 



new superposition rule to obtain the general solution for the corresponding 
Milne-Pinney equation is given. 

Another interesting question to be solved later on is why the obtained 
superposition rule for the Milne-Pinney equation and the already known one 
just depend on the variables X\ and %2 and not on the V\ and Vi variables as 
one would expect by the theory of Lie systems. This fact is not an isolated 
result and can be found also in other superposition rules for second order 
differential equations, see [13] . 
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